MPC=SS_MPC(PP,SS);
Temp_1  =   sum(reshape(MPC.QQ.*MPC.QW,SS.DistApp.fzo.UnitNum'.*[1,2]),2);
Temp_2  =   sum(reshape(MPC.QW,SS.DistApp.fzo.UnitNum'.*[1,2]),2);
AvgMPC  =   Temp_1./Temp_2;


%% Figure
% Graphics
N_col       =   1;
N_row       =   1;
fig_size    =   [N_col,N_row].*[1,0.75]*1/2.5; 
fig_gap     =   [0.02,0.02];
fig_VMargin =   [0.15,0.02];
fig_HMargin =   [0.18,0.02]; 

fig_SubPlot =   @(ii_vv)subtightplot(N_row,N_col,ii_vv,...
                                     fig_gap,fig_VMargin,fig_HMargin);
LEG         =   @(PH,LabelList)legend(PH,LabelList,'fontsize',8,'interpreter','latex','location','best');
XLabel      =   @(Label)xlabel(Label,'fontsize',8,'interpreter','latex');
YLabel      =   @(Label)ylabel(Label,'fontsize',8,'interpreter','latex');

AxisSetup   =   GR_Axis('TickFmt',struct('X','%g','Y','%.2g'),...
                        'FontSize',struct('X',8,'Y',8),...
                        'Exponent',struct('X',0,'Y',0),...
                        'Box','off');

Tick_X      =   @(X)GR_TickValue('Unit',4).Setup(X);
Tick_Y      =   @(Y)GR_TickValue('Symmetric',0).Setup(Y);


LP          =   GR_Line('Style','-','Color','Black','Width',3);
Fig         =   FigureSetup(['Average MPC by Wealth'],fig_size);
ax          =   fig_SubPlot(1);
LP.Plot(ax,SS.DistApp.fzo.UnitNode{1},AvgMPC);

XTickValue  =   [-1,0,1,5,10];%GR_TickValue('Unit',20).Setup(SS.DistApp.fzo.UnitNode{1});
YTickValue  =   [0,0.1,0.2,0.3,0.4,0.5];
ax.YTick    =   YTickValue; ax.YLim = [min(YTickValue),max(YTickValue)];
ax.XTick    =   XTickValue; ax.XLim = [min(XTickValue),max(XTickValue)];
AxisSetup.Setup(ax);
XLabel('Net Wealth (r.t. steady-state avg. qtr labor income)');
YLabel('MPC');
print('-depsc','-r1000',Fig,['TableGraphs\\Baseline\\MPC_Wealth']);
% close(Fig);
